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COMPARISON  OF  THE  ROUGH  SURFACE  REFLECTION  COEFFICIENT 
WITH  SPECULARLY  SCATTERED  ACOUSTIC  DATA 


INTRODUCTION 


Miller  and  Vegh  (1]  in  treating  reflection  from  the  rough  surface  of  the  sea  derived  a  one- 
parameter  family  of  curves  for  the  rough  surface  reflection  coefficient  or  roughness  factor  R  given  by 

R  (A’,e)  =  exp  [— 2e^T)^  (27r/?)^l  /q [26^ T)^(27r/f 1 


-t-  (1  —  e  ^  '^^exp  [— 4t) ^(27rj?) 

-  ^e^(l  —  €^)<I>|[-y,  1;  2;e^  -  {1-ng)^]  (1) 

where 


and 


g  =  (ctI\)  sin 


r,  =  IJ  ^(1  -*2)]-i/2 


Here  g  is  a  measure  of  the  effective  surface  roughness  or  simply  surface  roughness,  e(0  <  e  <  1)  is 
the  spectral  width  parameter,  <r  is  the  standard  deviation  of  the  water  surface  elevation,  is  the  grazing 
angle  for  specular  reflection,  A  is  the  wavelength  of  the  incident  radiation,  and  /o(-^)  is  the  modified 
Bessel  function  of  order  zero.  The  function  <t>|  (a,  jS;  y;  x,  >' )  is  a  confluent  hypergeometric  function 
in  two  variables  first  defined  in  1920  by  P.  Humbert  l2,  p.  58).  In  the  Appendix  we  derive  an  integral 
representation  for<l>i  that  may  be  used  for  numerical  computation. 

R  (g,  f),  given  by  Eq.  (1),  is  essentially  the  Fourier  transform  of  the  probability  density  D(y,i) 
for  surface  elevation  y  where 


DCv,6)  = 


27r^/' 


T)(T 


exp  ( 


(1  - 
TT^^^Tjrr 


exp  ( 


-V" 


4t)^ 


)  {cos  '  €  -f  €  (  1  —  6  A'  (26  ^ 


(2) 


Here  Aglx)  is  the  MacDonald  function  or  Bessel  function  of  imaginary  argument  of  order  zero. 
Kfja,  x)  is  an  incomplete  Lipschitz-Hankel  integral  of  A',)(x)  and  may  be  written  in  closed  form  either 
in  terms  of  incomplete  cylindrical  functions  131  or  in  various  ways  in  terms  of  Kampc  de  Feriet  func¬ 
tions  14,5];  e.g. 

,  r )  =  r  A'olr  )  A^(a  ,  z)  +  z^  K\(z)  /!()(«  ,  r  ) 
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where 


AAa,z)  =  F 


0:1;1 

2;0;0 


1/2;  1;  ,2 

I/2,3/2:  T 


I  0;2;1 
+  ^2:1;0 


-■  1,1; 

1,2;  3/2;  -;  T 


A^ia ,  z)  =  F 


0:1;1 

2:0;0 


1/2;  1;  ^2^2  ,2 

3/2,3/2:  -;  -;  T 


2; 

2:1;( 


2,2;  3/2;  -;  4  ’  4 


e),  given  by  Eq.  (2),  was  derived  in  Ref.  1  by  assuming  that  the  water  surface  could  be  described 
locally  by  sinusoids  with  uniform  phase  distribution  whose  amplitude  distribution  is  given  by  a  density 
function  derived  by  Rice  [6]  and  by  Cartwright  and  Longuet-Higgins  (7).  Figure  1  gives  graphs  for 
D(y,e),  for  various  values  of «. 


Fig.  1  -  Density  function />(>•,<)  lor  various  values 
of  the  spectral  width  parameter  c 


COMPARISON  OF  R  (g,€)  WITH  ACOUSTIC  DATA 

In  1980  DeSanto  (8,  p.  70,  Fig.  5|  compared  1)  with  acoustic  data  fnim  ('la>.  Mcciwin.  anti 
Wright  (9|.  Although  Rif;.  1)  was  first  derived  in  1974  IIO),  a  mathematicall\  tigorous  derisation  vsas 
not  obtained  until  1984  (11).  In  view  of  l;q.  (I),  it  now  appears  appropriate  to  compare  Ri,K.t)  with 
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the  aforementioned  data.  Whereas  1)  takes  into  account  only  the  standard  deviation,  ir.  oi  sur¬ 
face  elevation,  is  dependent  on  e  also  and  hence  on  the  moments  ol'  the  frcqucncv  energy 

spectrum  (.s )  of  the  surface  through  the  equations  (12.  p.  2461 

€  -  =  (m„  W4  —  m2  h' >>h)  '^'4 


Joe 

)  Js  (m,,=  it) 

0  " 

Figure  2  compares  1/3)  with  the  data  given  by  Clay  et  al.  in  Fig.  5  of  Ref.  9.  1/3) 

appears  to  be  in  better  agreement  with  this  data  than  the  multiple  scattering  theoretical  result  given  in 
Fig.  5  of  Ref.  8. 


Fig.  2  —  Comparison  of  (he  theoretical  curve  1/3) 
with  experimental  data 


CONCLUSION 

One  of  the  family  of  rough  surface  reflection  coefficients  agrees  with  acoustic  data  reasonably 
well;  at  least  as  well  as  the  curve  given  previously  by  the  multiple  scattering  model 
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Appendix 

INTEGRAL  REPRESENTATIONS  FOR  «1>|  [a.p,y\x,y] 


The  confluent  double  hypergeometric  function  <1>|  is  defined  by 


^  r  I  _  ^  +  nt  X"'  y 

<t>,  -y;  =  £ - Pi - 

/»,«  =  ()  'yhn*n 


m  y,n 

I  «  ' 


U-  I  <  I  ,  l.v  I  < 


The  definition  of<t>|  given  in  Erdelyi  et  al.  (Al,  p.  225]  and  Gradshteyn  et  al.  [A2,  9.261,  Eq.  1]  is 
incorrect. 


By  using  Ref.  A3,  p.  266 


(« )„ 


r  (>) 


r>+«*-  I 


(y)^  r(a)r  (y  -  a)  0 
with  the  definition  of  <t>|  given  above  we  obtain 


(1  -  Rey  >  Re«  >  0 


<t‘  I  (o ,  /i ;  y ;  .v ,  ,y  1  = 


IhL 


Z  Z  /  '  t'”*”*"-'  (1  -  (fi) 


r(a)r  (y  -  a) 


m-0  »i"»0 


m\n\ 


dt 


Now  interchanging  the  integral  sign  and  double  sum  and  noting  that 


(fx)" 


n-0 


m-0 


W  ! 


we  obtain  for  Rey  >  Rea  >0,  U|  <  1,  |y|  <  oo 

d>  1  [a ,  /3 ;  y ;  .V ,  w  I  =  — — ^ - -  f  e”  { 1  -  .vt  l  ( 1  -  r  " '  t"'  '  df 

I  (a  )  I  (y  —  a )  o 


In  particular. 


<b,  (3/2,  1;  2.x,>]  =  —  f  ' 

77  (I 


(,0,1/2 


dr 


n  ''  n  (1  _  x7)  (1  - 

Now  making  the  transformation  t  =  sin*  H  and  replacing  v  by  and  >  by  -  t’r’  we  obtain 
<l>ii 


,13/2,  1;  2;t^  =  -  f" 

'  77  " 


4  Sin‘«  (' 

-  uH 


I  -  « ’  sin’ W 


(Al) 


for  real  t,  >  the  integrand  here  is  nonnegativc  on  the  closed  interval  |0.  7r/2|  and  haN  no  singularities 
for  0  ^  <  I;  the  integral  in  Eq.  (Al)  is  therefore  suitable  for  numerical  quadrature  and  4>|  may 

thereby  be  computed. 
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It  may  also  be  shown  [I,  Eq.  15]  that 


4>,  13/2,  1;  2:e^ 


=  2,2]  = 


2(1  -  e^) 


2)1/2 


2  /(^-''y„(2>7)  erf 


(1  -  *2)1/2 


dt 


from  which  it  follows  that 

lim  *2  (1  _  ^2)  ^  [3/2^  I;  2;  *2^  — ^2^,3]  ^  Q 

«  —  1 


(A2) 


Hence  Eq.  (1)  is  valid  in  the  limit  for*  =  1. 
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